Damped collective motion of isolated many body systems within a variational 

approach to functional integrals 
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Two improvements with respect to previous formulations are presented for the calculation of the 
partition function Z of small, isolated and interacting many body systems. By including anhar- 
monicities and employing a variational approach quantum effects can be treated even at very low 
temperatures. A method is proposed of how to include collisional damping. Finally, our approach 
is applied to the calculation of the decay rate of metastable systems. 

PACS numbers: 24.10.Cn, 05.30.-d, 05.60.Gg, 82.20.Xr 



Small many body systems like atomic nuclei and metal 
clusters may undergo self-sustained collective motion. In 
the former case this has long been established experimen- 
tally for situations where the nucleons may be assumed 
to move in a time-dependent state with no additional ex- 
citation energy j]|. In the last few years one is also able 
to measure properties of single, isolated metallic grains 
y . Although such systems consist of (many) particles of 
identical or similar nature collective motion only involves 
one or a few time dependent parameters. The latter must 
be introduced so that they obey certain basic relations 
for the intrinsic (particle) degrees of freedom. One, for- 
mally promising way is to apply the functional integral 
technique and to introduce collective degrees of freedom 
via a Hubbard- Stratonovich transform (HST). This is es- 
pecially useful when the nature of the collective motion 
can be guessed and attributed to one generator F. Fol- 
lowing 1] one may then introduce a separable two-body 
interaction to obtain for the total Hamiltonian 



H = Ha 
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FF, 
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in which the coupling constant k is negative for iso- 
scalar modes. Of course, this might be understood as 
one prominent term in an expansion of a general two- 
body interaction V^"^^ into a complete series of separable 
terms like H ^f + T>(2) = Hq + ^ TJ^=i k,, F'^F'^. For 
the sake of simplicity we will restrict our discussion to the 
model case ^ , but treatments of more terms are feasible 
and have already been undertaken for applications which 
are simpler than those we have in mind (see e.g. j^])- 
For zero thermal excitation the Hq might be assumed to 
simply represent the dynamics of independent particles. 
However, in case the constituents themselves are excited 
this may no longer be true as then incoherent, residual 

" (2) 

interactions V^-cs may play a role, with an ever increas- 
ing role the larger the excitation. In the nuclear case it 
is known that single particle excitations acquire a finite 
width if they are away from the Fermi level by only 5 — 8 
MeV y . This effect is further enhanced if the system of 
particles is thermally excited. In the present Letter we 
want to show how these effects may be accounted for in 



the functional integral approach. Moreover, in contrast 
to previous formulations [3|, |5|, |y{, we also demonstrate 
the usefulness of generalizing the variational formulation 
of |l2, iS] to many body systems. In this way one may 
overcome the divergence known to exist for small tem- 
peratures in the case where the collective motion passes 
over barrier regions. 

Our basic theoretical tool will be the partition func- 
tion Z which at given temperature T = l/(fcB/5) may be 
expressed as |3| 
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The J-^^^{P,qo) is a free energy which depends on the 
collective variable go, treated here in the static limit, 
which in this context in the literature is referred to 
as "Static Path Approximation" (SPA) '5'|. It may 
be written as T^'^^{P,qo) = '7o/(2|fc|) - Tin z{l3, go), 
where z(/3, go) is the grand canonical partition func- 
tion calculated for the static one body hamiltonian 
TinSTbo] = ^0 + QoF- Actually, the collective variable 
is "time-dependent" and is introduced via the HST. Us- 
ing this manipulation the two body interaction kFF/2 
disappears. The dependence on the imaginary time 
T is treated through the Fourier series ^(t) — qa + 
l^r^o 9'" exp(ii^r''"), with the Matsubara frequencies Vr = 
2TTr/{hf3) — 2TTrT/h. Genuine quantum effects in collec- 
tive motion are hidden in the factor 



C(/3,9o) 



/ V'q exp[-SE(/3,go)/^]- 
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The key stone for improvements over the SPA is given 
by the Euclidean action se- In the so called "Perturbed 
SPA" (PSPA) - also known as "SPA+RPA" or "Corre- 
lated SPA " (CSPA) - this Se is expanded to second order 
in the qr- In this way quantum effects are treated at the 
level of a local RPA, see e.g. ^^, ^J. In an extension of 
this approximation we shall want to make use of terms up 
to fourth order and write the action as (see also |l(lll ij) 






Prst qrQsqt + ^ cTrstu qrqsqtqu I + 0{ql 

r,s.t,u^O 
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The interesting point is that the coefficients A, p and a can be expressed by the one body Green functions associated 
with the Hamiltonian TiHSTi^o]; which at first may be assumed to represent simply independent particles for which 
the Green function is gj^ (z) = (z — ek)^^ ■ In such a case the coefficient (Trstu of fourth order consists of terms like 
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to give just one example of what these quantities look 
like. 

At the level of PSPA it is useful to connect the coef- 
ficient A to the response function_x(i^) defined through 
5{F)u, = -X{^)5q{^)- One gets [^ 
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The Xr serves as the stiffness in g^-direction: 
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The unperturbed intrinsic excitations hujki — ^k ~ ^i are 
to be calculated from the eigenvalues of HhstI^o]- The 
137^, on the other hand, are to be found from the secular 
equation 1 + fc xC^a*) = ^^'^ represent the local RPA 
frequencies. For an unstable collective mode /i the vu'f^ is 
negative such that Ai may become negative as well. This 
happens for temperatures below the so called crossover 
temperature Tq. There the dynamical fluctuations in gi- 
direction become too large for the harmonic approxi- 
mation to be justified. Formally, this shows up when 
in (O Gaussian integrals become divergent [3, U^ uM- 
The name crossover temperature is borrowed from work 
on "dissipative tunneling" within the Caldeira-Leggett 
model (CLM), where a similar breakdown is observed, 
see e.g. [ijl. From this work it is known that in the 
crossover region T w Tq the situation can be cured by 
treating the anharmonicities in (ji -direction explicitly up 
to fourth order. In this way it is possible to remove the 
divergence of the decay rate at Tq in the CLM [13 . In 
|10| this technique has been applied to many body sys- 
tems to overcome the shortcomings of the PSPA. Using 
this "Extended PSPA" (ePSPA) the partition function 
can be evaluated down to T = Tq/2 where the fluctu- 
ations in (72-direction become large, too. At very low 



temperatures the harmonic treatment of the fluctuations 
is no longer justified in any direction q^- This feature 
prevents analytic treatments of higher modes g^- 

Such problems may largelybe abolished by applying a 
variational procedure. In |2|, |a| the latter has been devel- 
oped for the system of a particle moving in a one dimen- 
sional potential. It allows one to evaluate the quantum Z 
at arbitrary temperatures. Using the formalism of coher- 
ent states, these ideas have been applied to many body 
systems in 'l4J. Here we want to make use of the expan- 
sion Q). Details of this novel method will be given in a 
forthcoming paper. The main idea consists in rewriting 
® as 

C(/3,go) = Jv'q exp[-4V?i] exp[-(sE - 4,°)/^] 
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The reference action s^ is introduced to specify an aver- 
aging procedure for which one may exploit the inequality 



(exph(sE-4°)A])a >exph(sE 
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for an optimization of the expression on the right hand 
side, which is easier to evaluate than the average in ||S}. 
The reference action should be chosen such that the nor- 



JV'q exp[-s1^/h] in © can 



malization factor Cq" = 
be evaluated exactly. A reasonable choice can be con- 
structed from the PSPA action Q by replacing the RPA 
frequencies nj^ by variational parameters fi^ 
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Evidently, when calculating the integrals for (se — Sfj") 
all terms disappear which are odd in g^- Hence, in the 
truncated expansion Q only terms of second and fourth 



order survive. Using the abbreviation II,, = Hul^r 
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FIG. 1: Free energy of the many body system for the LMGM 
with effective coupling constant « = 1.313 in various approx- 
imations. The arrow points to the inverse cross-over temper- 
ature /3o = l/To where the PSPA breaks down. 



Collecting all contributions the dynamical corrections 
read 
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which cannot be larger than In^. The advantage of this 
novel method over the PSPA and the ePSPA is that it can 
be applied at any T . The reason is that in contrast to the 
case of the PSPA the factors i/^ + i^'i, which contribute 
to the stiffness in g^-direction in 1)10(1 , stay positive for all 
n^. Consequently all Gaussian integrals in 1)13(1 converge. 
To test the results we apply the exactly solvable 
Lipkin-Meshkov-Ghck model (LMGM) 15], as used be- 
fore in 110 for the the SPA and the PSPA. In Fig.[Tlthe 
free energy ^ = — T In Z associated to the total Hamilto- 
nian (^ is compared with the approximations mentioned 
before. For the calculation the same set of parameters has 
been taken as in [iflj . The classical SPA gives reasonable 
results only at small (3 (high T). Inclusion of quantum 
effects at the level of local RPA within the PSPA delivers 
considerable improvement at not too large /3. However, 
even before its breakdown ai P — f3o it becomes unre- 
liable. The ePSPA behaves completely smoothly in the 
crossover region but breaks down at /3 = 2/3o. Compared 
to these results, the improvement found for the varia- 
tional approach is very striking. Notice please that it is 



free of discontinuities and the relative error is only of the 
order of 1%, even at very large /3. 

So far the one body Green's functions g^ (z) have been 
calculated within the picture of independent particle mo- 
tion. As mentioned before, within this model the residual 

' (2) 
two body interaction T^es is neglected. This interaction 

describes the incoherent scattering of particles and holes 
and couples Iplh states to more complicated ones. This 
mechanism may be understood as the origin of the damp- 
ing of collective motion, see e.g. [l^. In the following we 
will account for such couplings by dressing the one body 
Green functions with self-energies E = S' — ir/2 meaning 
that g^. (z) is replaced by 



gPi^) 



(z-efe-EUz)+ir,(z)/2)- 
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The dependence of the width Tk{w) on frequency and its 
variation with temperature is parameterized by the form 



2t^2 



Tk{uj) = 
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suggested in IT]. For zero temperature it is in good 
agreement with empirical data for the widths of proton 
and neutron states [j| . Finally, some comments are in or- 
der concerning our approximation in handling the impact 

" (2) 
of 14os • What is definitely implied in evaluating many- 
body Green functions is the application of a factorization 
assumption. However, it should be noticed that this is 
done for excitations of the intrinsic dynamics and not for 
the collective modes. One may therefore argue that co- 
herent effects are small, in particular at larger thermal 
excitations. After all the SPA is meant to represent the 
high temperature limit. 

Having established the connection between the PSPA 
and the theory of linear response one may take over the 
method developed there (see [la|) to extract transport 
coefficients for collective dynamics. For slow collective 
motion, as given for nuclear fission, it suffices to con- 
centrate on the lowest mode of the collective response 
function. In that frequency regime, the latter may then 
be approximated by the response function of a damped 
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FIG. 2: A comparison of different approaches to the decay 
rate for a particle in a one-dimensional cubic potential with 
barrier, for zero (77 — 0) and critical (r; = 1) damping. The 
arrows point to the /3o — I/Tq where the PSPA breaks down. 



harmonic oscillator, which means the replacement: 



Xcoll(w) 



x(^) 



1 + kxH 



Xosc(w) (16) 



At any go the Xosc(w) is uniquely related to the coef- 
ficients of the local motion such as the inertia, friction 
and stiffness. It should be mentioned that exactly at 
this point an important difference to the CLM shows up. 
The latter only allows a microscopic interpretation of dis- 
sipation but not for the collective potential nor for the 
inertia. Conversely, our approach allows us to determine 
all coefficients in a self-consistent fashion, and all of them 
are influenced by the microscopic damping mechanism. 
It may be noted that the replacement of 5^ (z) by the 

gf, (z) also influences the anharmonic terms of the Eu- 
clidean action, see (0. With respect to the variational 
approach the reduction (|16|l leads to a considerable sim- 
plification: Only one variational parameter fl has to be 
dealt with in ^ to (jTSJ. 

Let us finally turn to the decay rate of meta-stable 
states. As shown by J. S. Langer [13 this rate can be 
related to the imaginary part of the free energy. Indeed, 
the CLM makes extensive use of this feature [13. We 
will take this procedure over to the case of the interact- 
ing many body system. In order to obtain Im J- from the 
2 in (|21) the integration contour has to be deformed into 
the complex qo-plane. This has to be done at the barrier 
of J^^^^ in the sense of a steepest descent approximation. 
Details of the resulting rate formula for interacting many 
body problems will be given in a forthcoming paper. The 
quantum corrections /qm = Cb/Ca > 1 to the purely clas- 
sical rate formula of the SPA can be calculated from the 
dynamical factor C of (O evaluated at the barrier b and 
at the minimum a of JF^ within the approximations 
discussed above. 

To compare the various approaches to Z with respect 
to the evaluation of the quantum decay rate we take the 



case of a particle moving in a one dimensional cubic po- 
tential with barrier |l3l |. For such a situation our varia- 
tional approach reduces to the Feynman-Kleinert method 
(FKV) 8]. Results are shown in Fig. 13 The classical rate 
(SPA) just displays the exponential behavior of the Ar- 
rhenius law. Quantum corrections are known to increase 
the decay rate at larger /3 (low T). The variational ap- 
proach smoothly matches the classical result to the low 
T result obtained by Grabert et al. with the dynamical 
"bounce solution" jl3|. In the region /? < /So the PSPA 
would show a similar behavior with the exception of the 
crossover region itself. Different to the PSPA the vari- 
ational approach gives good results also for /3 « /3o but 
also fails in the regime (3 ^ (3^. This behavior may per- 
haps be understood as follows. Our variational approach 
in its present version starts from a static approximation 
and only includes small quantum corrections. For very 
large (3 (low T) this cannot be enough as one ought to 
start from genuine time- dependent paths ^(t), like those 
associated to the "bounce solution" in the CLM. For the 
many body problem with dynamics beyond independent 
particle motion, however, this is still an open problem. 

The authors would like to thank J. Ankerhold and F. 
Ivanyuk for helpful discussions. 
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